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ABSTRACT 

A new class of solutions for Einstein’s field equations representing a static spher¬ 
ically symmetric anisotropic distribution of matter is obtained on the background of 
pseudo-spheroidal spacetime. We have prescribed the bounds of the model parameters 
k and po on the basis of the elementary criteria for physical acceptability, viz., regu¬ 
larity, stability and energy conditions. By taking the values of model parameters from 
the prescribed bounds, we have shown that our model is compatible with the observa¬ 
tional data of a wide variety of compact stars like 4U 1820-30, PSR J1903+327, 4U 
1608-52, Vela X-l, PSR J1614-2230, SMC X-4 and Cen X-3. 

Subject headings: General relativity; Exact solutions; Relativistic compact stars. 
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1. Introduction 

The study of compact objects in agreement with observational data has received wide attention 
among researchers. A number of superdense star models, compatible with observational data, 
have appeared in literature in the recent past ( Murad (2013a), Murad and Saba (2013b,c, 
2014a), Maurya et al. (2015) & Sharma and Ratanpal (2013)). If spacetime admitting compact 
star models possess a definite three-space geometry, then it is a mathematically interesting 
problem also. The spheroidal spacetimes studied by Vaidya and Tikekar (1982), Tikekar 
(1990) and the paraboloidal spacetime studied by Finch and Skea (1989), Tikekar and Jotania 
(2005), Sharma and Ratanpal (2013) are examples of spacetimes with definite 3-space geometry. 
The superdense star models developed by Tikekar and Thomas (1998) has pseudo-spheroidal 
geometry. A number of researchers used this spacetime for developing physically viable models 
of compact stars under different assumptions on the physical content. 


Theoretical investigations of Ruderman (1972) and Canuto (1974) suggest that matter may 
not isotropic in ultra high density regime. After the publication of the work of Bowers and Liang 
(1974), there has been a large number of models devoted to the study of anisotropic distribution of 
matter. Maharaj and Maartens (1989) developed an anisotropic model with uniform density and 
Gokhroo and Mehra (1972) gave a more realistic anisotropic model with non-uniform density. 
Tikekar (1999); Tikekar and Thomas (2005), Thomas et al. (2005) developed superdense 
anisotropic distributions on pseudo-spheroidal spacetimes. Thomas and Ratanpal (2007) studied 
non-adiabatic gravitational collapse of anisotropic distribution of matter accompanied by radial 
heat flux. Dev and Gleiser (2002, 2003, 2004) have studied the impact of anisotropy on the 
stability of stellar configuration. Anisotropic distributions of matter incorporating linear equation 
of state have been studied by Sharma and Maharaj (2007), & Thirukkanesh and Maharaj 
(2008). Komathiraj and Maharaj (2007) have studied charged distribution using linear equation 
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of state. Sunzu et al. (2014) studied charged anisotropic quark stars using linear equation of 
state. Anisotropic distributions of matter incorporating quadratic equation of state have been 
given by Feroze and Siddiqui (2011) & Maharaj and Takisa (2012). Varela et al. (2010) used 
linear and non-linear equations of state for describing charged anisotropic distributions of matter. 
Paul et al. (2011) have shown, in the MIT bag model of quark stars, that anisotropy can affect the 
bag constant. Polytropic equations of state has been used by Thirukkanesh and Ragel (2012) & 
Maharaj and Takisa (2013b). Malaver (2013a,b, 2014) and Thirukkanesh and Ragel (2014) have 
used modified Van der Waals equation of state for describing anisotropic charged compact stars. 


Recently Pandya et al. (2015) have developed anisotropic models of compact stars compati¬ 
ble with observational data by generalizing Finch and Skea (1989) ansatz. The anisotropic stellar 
model given by Sharma and Ratanpal (2013) is a subclass of the model of Pandya et al. (2015). 
This model accommodates the observational data of a variety of compact objects recently studied 
by researchers. In the present article, we have obtained a new class of anisotropic stellar model of 
compact objects on the background of pseudo-spheroidal spacetimes. The physical parameter po 
and geometric parameter k appearing in the model are restricted as a result of various physical 
acceptability conditions imposed on the model. Another geometric parameters R of the model 
plays the role of the radius of the spherical distribution of matter. It is found that our model 
yields values of different physical quantities that are in good agreement with the most recently 
available observational data of compact objects ( Gangopadhyay et al. (2013)) like 4U 1820-30, 
PSR J1903+327, 4U 1608-52, Vela X-l, PSR J1614-2230, SMC X-4 and Cen X-3. 


We have organized the paper as follows: In section 2, we have solved the field equations and 
obtained restrictions on the model parameters using various physical requirements and energy 
conditions. The bounds for the model parameters k and po are obtained in section 3. In section 


-5- 


4, we have shown that our model is compatible with recent observational data of a number of 
compact objects ( Gangopadhyay et al. (2013)). The main results obtained in the present work is 
discussed in section 5. 


2. Spacetime metric 

A three-pseudo spheroid immersed in four-dimensional Euclidean space has the Cartesian 
equation 


u 2 x 2 +y 2 +z 2 _ 
b 2 ~R 2 

The sections u = constant are spheres of real or imaginary radius according as u 2 > b 2 or 
u 2 < b 2 , while the sections x = const , y — const, and z — const are respectively, hyperboloids of 
two sheets. 

On taking the parametrization 

x — Rsinhksindcosty 
y — RsinhXsinQsinO 
Z — RsinhXcosO 

u — bcoshX (1) 

the Euclidean metric 

do 2 — dx 2 + dy 2 + dz 2 +du 2 

takes the form 

1 4- k -t 

do 2 — - ^dr 2 + r 2 dd 2 + r 2 sin 2 Gd(j) 2 


( 2 ) 
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/ 2 

where k= 1 + fi and r = RsinhX. The metric (2) is regular for all points with k > 1 and call 
pseudo-spheroidal metric (Tikekar and Thomas (1998)). 

We take the interior metric describing the anisotropic matter distribution in the form 



( 3 ) 


where, k , R are geometric parameters and k > 1. This spacetime, generally known as pseudo- 
spheroidal spacetime, has been studied by many researchers (Tikekar and Thomas (1998); 
Tikekar (1999); Tikekar and Thomas (2005); Thomas et al. (2005); Thomas and Ratanpal (2007); 
Paul et al. (2011); Chattopadhyay and Paul (2010); Chattopadhyay et al. (2012)). 

Following Maharaj and Maartens (1989), we write the energy-momentum tensor for anisotropic 
matter distribution in the form 


T ij = (P+p) u i u j ~ PSij + Xij, 


( 4 ) 


where, p,p and u l denote the proper density, fluid pressure and unit four-velocity of the fluid, 
respectively. 

The anisotropic stress-tensor 7l t j is given by 


tt'ij — y/3S CjCj ( UjUj gij ) , 


( 5 ) 


where, C l — (0, — e _/l / 2 ,0,0) is a radial vector and S — S(r ) denotes the magnitude of the 
anisotropic stress. 

The non-vanishing components of the energy-momentum tensor are given by 



( 6 ) 


Hence the radial and transverse pressures are given by 


Pjl 


Pr 



( 7 ) 


( 8 ) 
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Then the magnitude of the anisotropic stress has the form 


5 = 


Pr~P± 

V3 ' 


(9) 


The physical and geometric variables, related through Einstein’s field equations 




( 10 ) 


are to be determined from the following set of three equations: 


8np = 


8k p r — 


l-e~ x e~ x k' 


+ ■ 


e A - 1 e A v 


r 

- 2 ,,/ 


8np± = e 


r 

v" v /2 v'A' V-k' 
2 + 4 4 + 2r 


(ID 

( 12 ) 

(13) 


where a prime denotes a differentiation with respect to r. The equations (11) - (13) can be 
couched in the form 


2m \ , 

l ~T) v 

~(SnV3S) 


= 1 - 


2m 


2m 

8np r r H— 


= (8^p + 8^p r )v / + 2(8^p' r ), 


where 

r 

m(r ) = 47T j u 2 p(u)du. 

0 

The energy-density p and the mass in within the radius r have expressions 


(14) 

(15) 

(16) 

(17) 


8k p 


k I 


3 + 


/?“ / 2 ^ ^ 


(18) 


m(r) = 


* (*-!)£ 
2 1+^4 

^ w 


(19) 
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It can be easily obtained from equation (18) that 


8 Tip' = 


2 k{k — 1 )r 


5 + k 


RV 


R 4 


l + k W 


< 0 , 


( 20 ) 


indicating that the density decreases radially outward. 

In order to obtain the metric potential v, we assume an expression for p, in equation (15), in the 
form 


87 ip r — 


P 0 


1 


r 2 

w 


1 + k 


R 2 


l + k W 


( 21 ) 


The radial pressure in the present form vanishes at r — R and takes the value ^4 at the centre r = 0. 
It is non-negative for all values of r in the range 0 < r < R. Further, on differentiating equation 
(21) with respect to r, we get 


8 Kp r — 


Ap 0 r [k+jp 

^ 4 (l+k| 


< 0 , 


( 22 ) 


indicating that the pressure p, decreases radially outward. Since p r (r — R) — 0, the geometric 
parameter R takes the role of the boundary radius of the distribution. With this choice of p r , 
equation (15) can be integrated to obtain v in the form 


e v =A \l+k 


pp(k+l) 
" 2 \ 2k 2 


R 2 

exp 


k -1 

.2 \ ~ 


1 + 


R 2 


-Po 
2 k 2 


l+k w 


(23) 


where A is a constant of integration. 

Therefore, the spacetime metric takes the explicit form 


ds — A ( 1 T 


PQ{k+l) 

2k 2 


- W 

r 2 \ 2 

1 + F 


xexp 


~P 0 
2k 2 


1 +k 


R 2 


dt 2 - dr 2 — r 2 d0 2 — r 2 sin 2 0d(j) 2 . 


1 + 


(24) 


W 
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The constant of integration A can be obtained by matching the interior spacetime metric (3) with 
the Schwarzschild exterior metric 

1 — —^ dr 2 — r 2 d0 2 — r 2 sin 2 0d(j) 2 (25) 

across the boundary r — R. This gives 



R 


2M{k +\) 
(*- 1 ) 


(26) 


and 


A 


2 

k +I 



(*+i)po 



(27) 


The expression for anisotropy is now readily available by substituting for p r , p' r and v' in the 
equation (16). 


8 7tV3S 


^ 2p 0 


B(r)C(r) 

\ 

U ( 

! + ^) 3 

Rl ( 1+k B / 



(28) 


where 


B(r) = 


Po\ l ~W) k -1 

+ ■ 


4 1 + S 4 1 + F 


C{r) = 


(k~ 1 ) ( 3+kjp ) +Po ( 1 - 


R 4 


(29) 

(30) 


It is easy to see that S vanishes at origin r — 0, which is a desired requirement for anisotropic 
distributions ( Murad (2013a), Murad and Saba (2013b,c, 2014a) & Bowers and Liang (1974)). 
The expression for transverse pressure 


87 ip±_ — 87 ip r — 8 7tV3S 


(31) 


can be obtained using equations (21) and (28). 

Moreover, the condition p± > 0 will lead to the following inequality at r — R 

Po < J^( k ~ !) 2 (^ + 3), 


(32) 
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whereas at r = 0 , the condition is evidently satisfied. 

dp,- nnr | dp± 
dp ana dp 


The expressions for ^ and c ^~ are given by 


dp r 2po 




dp k(k — 1 ) 


5 + k£ 


dpi. dp r -dS 


dp dp 


dp 


The conditions 0 < ^ < 1 and 0 < < 1 at r — 0, respectively, give the inequalities 


dp± 


dp 


dp 


0 <Po< 


and 


0<po< 


k(k-l)(k + 5) 


2(k+l) 

Similarly the above conditions at r = R, respectively, give 

2(3 k+ 1) - 33k 2 + 30k + 1 < p 0 < 2(3 k+ 1) - ^12>k 2 + 50k+1, 


(33) 


(34) 


(35) 


(36) 


(37) 


and 


The adiabatic index 


has the explicit expression 


0 <po < 


r = 


-k 4 +12k 3 + 78k 2 -92k+ 3 
8k 2 — 24k + 80 ' 

r _ P +Pr dp r 
Pr dp 

2 (k+ jp) C(r) 


k(k-l)[5 + k£) l-£ 


(38) 


(39) 


(40) 


The necessary condition for the model to represent a relativistic star is that r > 3 throughout the 
star, r > 4 at r — 0 impose a condition on po, viz., 


k -1 


Po > 


3 


(41) 
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The strong energy condition p — p r — 2pj_ > 0 at r — 0 and r = R, respectively, give the following 
two inequalities 

Po < k - 1 (42) 


and 

^ (* + 3)(*-l)(*-5) 

Po - - 16 - 

In order to obtain a valid range for the parameters po and k, we have to consider the 
inequalities (35) - (38) and (40) - (43) simultaneously. 


(43) 


3. Bounds for Model Parameters 

The pseudo-spheroidal space-time model developed for anisotropic matter distribution contains a 
physical parameter po related to the central pressure and two geometric parameters, viz., R and k. 
Since p r [r — R) — 0, the free parameter R represents the radius of the distribution. The bounds for 
the other two parameters po and k are to be determined by the following requirements a physically 
acceptable model is expected to satisfy in its region of validity, 0 < r < R. 


1 ) p(r) > 0, p r (r)> 0, p±(r)> 0; 


2 ) p(r) -p r (r) -2p ± (r) > 0 ; 


3 ) 


Mrl K o 

dr ^ 


¥<« 


°<t 


<i. 
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5) The adiabatic index T(r) > 


The conditions p(r) > 0,p r (r) > 0, < 0, - < 0 are automatically satisfied by equations 

(18), (21), (20), (22). 

We have displayed in Table 1 the bounds on po in terms of the parameter k at the centre and on 
the boundary. 


Table 1: Bounds for po. 


Physical requirements at r = 0 


at r — R 


P Pr 2 p ± >0 p 0 <k-l 


> (k+3){k- l)(A->) 

PO ^ 16 


0<^<1 0<p 0 < §(fc—1) 2(3 k+ 1) - ^(33k 2 + 30k+ 1) < p 0 

< 2(3 k+ 1) - ^(I3k 2 + 50k+l) 




0<po< 


k(k-l)(k+ 5) 
2{k+l) 


0<po< 


■k A : 12A- 3 +78/c 2 —92A.--I-3 

8 & 2 — 24&+80 


r(r) > f 


PO > 


/t-l 

3 


Automatically satisfied 


We have displayed the numerical values of the lower and upper bounds of po for different values 
of k > 1 in Table 2. We have considered the maximum of all lower limits of po and minimum 
of all its upper limits. The admissible values of k are those for which minimum of upper limit 
minus maximum of lower limit is positive. This condition restricts the values of k in the range 
(2.05,5.69). It is further observed that for 2.05 < k < 3.47, 3.47 <k< 5.24, and 5.24 <k< 
5.69, po satisfies, respectively, the inequalities < po < 11 , < po < k— 1 and 
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^(k + 3)(k— l)(fc — 5) < po < k— 1. The shaded region in Figure 1 gives the permissible values 
of k and pq. Any values of k and po outside this region may violate one or other of the physical 
requirements of the model. 



k 


Fig. 1.— Permissible values of po and k 
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Table 2: Permissible values of k and pq. 


Here, LPl = ZJ°2 — (*+ 3 K t ~ 1 )(*-5) ; LP3 = 2(3it+ 1) - y/(33k 2 + 30k+ 1), 


UPl = 

UP4~ 

. (K+3)(k-l) 2 

16 ’ 
k(k—l)(k+5) 

UP2 = 

UP5 = 

k- 1, 

2(3Jfc + 

UP3 = 

5(H) 

2 ’ 


UP6 = 

-U+12£ 3 +78fc 2 - 

-92k+3 


l)-^(l3k 2 + 50k+l) : 

' 2(k+l) ’ 

8 k 2 

—24£+80 

k 

Lower Limit for po 

Max 


Upper Limit for po 


Min 

Min - Max 


LPl 

LPl 

LP3 


UPl 

UPl 

UP3 

UP4 

UP5 

UP6 

2 

0.33 

-2.50 

0.11 

0.33 

0.31 

1 

2.5 

2.33 

1.63 

3.30 

0.31 

-0.02 

2.05 

0.35 

-2.61 

0.12 

0.35 

0.35 

1.05 

2.625 

2.49 

1.72 

3.54 

0.35 

0.00 

2.1 

0.37 

-2.71 

0.12 

0.37 

0.39 

1.1 

2.75 

2.65 

1.82 

3.78 

0.39 

0.02 

2.4 

0.47 

-3.28 

0.18 

0.47 

0.66 

1.4 

3.5 

3.66 

2.40 

5.32 

0.66 

0.19 

2.8 

0.60 

-3.83 

0.26 

0.60 

1.17 

1.8 

4.5 

5.17 

3.21 

7.40 

1.17 

0.57 

3 

0.67 

-4.00 

0.30 

0.67 

1.50 

2 

5 

6.00 

3.63 

8.40 

1.50 

0.83 

3.1 

0.70 

-4.06 

0.32 

0.70 

1.68 

2.1 

5.25 

6.43 

3.84 

8.88 

1.68 

0.98 

3.2 

0.73 

-4.09 

0.35 

0.73 

1.88 

2.2 

5.5 

6.87 

4.05 

9.35 

1.88 

1.14 

3.4 

0.80 

-4.10 

0.39 

0.80 

2.30 

2.4 

6 

7.79 

4.48 

10.23 

2.30 

1.50 

3.47 

0.82 

-4.08 

0.40 

0.82 

2.47 

2.47 

6.175 

8.12 

4.63 

10.53 

2.47 

1.64 

3.8 

0.93 

-3.81 

0.48 

0.93 

3.33 

2.8 

7 

9.75 

5.34 

11.79 

2.80 

1.87 

4 

1.00 

-3.50 

0.52 

1.00 

3.94 

3 

7.5 

10.80 

5.78 

12.46 

3.00 

2.00 

4.2 

1.07 

-3.07 

0.57 

1.07 

4.61 

3.2 

8 

11.89 

6.22 

13.05 

3.20 

2.13 

4.4 

1.13 

-2.52 

0.62 

1.13 

5.35 

3.4 

8.5 

13.02 

6.66 

13.58 

3.40 

2.27 

4.8 

1.27 

-0.99 

0.71 

1.27 

7.04 

3.8 

9.5 

15.41 

7.55 

14.45 

3.80 

2.53 

5 

1.33 

0.00 

0.76 

1.33 

8.00 

4 

10 

16.67 

8.00 

14.80 

4.00 

2.67 

5.2 

1.40 

1.15 

0.81 

1.40 

9.04 

4.2 

10.5 

17.97 

8.45 

15.10 

4.20 

2.80 

5.24 

1.41 

1.40 

0.82 

1.41 

9.26 

4.24 

10.6 

18.23 

8.54 

15.15 

4.24 

2.83 

5.4 

1.47 

2.46 

0.85 

2.46 

10.16 

4.4 

11 

19.31 

8.90 

15.35 

4.40 

1.94 

5.67 

1.56 

4.52 

0.92 

4.52 

11.82 

4.67 

11.675 

21.18 

9.52 

15.63 

4.67 

0.15 

5.69 

1.56 

4.69 

0.92 

4.69 

11.95 

4.69 

11.725 

21.32 

9.56 

15.64 

4.69 

0.00 

5.71 

1.57 

4.85 

0.93 

4.85 

12.08 

4.71 

11.775 

21.46 

9.61 

15.66 

4.71 

-0.14 
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4. Compact Star Models 

In order to validate the model, we examine our model with observational data. We have 
considered the pulsar 4U 1820-30 whose estimated mass and radius are 1.58M 0 and 9.1 km 
( Gangopadhyay et al. (2013)). If we set these values for mass and radius then from equation 
(26) we obtain the value of k — 3.1 which is well inside the valid range for k. Similarly 
assuming masses of some well studied compact stars like PSR J1903+327, 4U 1608-52, Vela 
X-l, PSR J1614-2230, SMC X-4 and Cen X-3, we have obtained the same radius calculated by 
Gangopadhyay et al. (2013) for values of k in the valid range. The values of mass, radius, k and 
other relevant quantities like central density p c , density at the boundary p«. the compactification 
parameter u and at the centre for po = 1.08. are shown in Table 3. 


Table 3: Estimation of physical values based on observational data for po — 1.08. 


STAR 

k 

M 

R 

Pc 

Pr 

«(=?) 

O 

II 

!<. 

"72 



(M q ) 

(Km) 

(MeV fm~ 3 ) 

(MeV fm~ 3 ) 



4U 1820-30 

3.100 

1.58 

9.1 

2290.97 

277.12 

0.256 

0.206 

PSR J1903+327 

3.176 

1.667 

9.438 

2129.82 

257.62 

0.261 

0.199 

4U 1608-52 

3.458 

1.74 

9.31 

2188.78 

267.75 

0.276 

0.176 

Vela X-l 

3.407 

1.77 

9.56 

2075.80 

251.08 

0.273 

0.179 

PSR J1614-2230 

3.997 

1.97 

9.69 

2020.48 

244.39 

0.300 

0.144 

SMC X-4 

2.514 

1.29 

8.831 

2432.67 

294.25 

0.215 

0.285 

Cen X-3 

2.838 

1.49 

9.178 

2252.20 

272.42 

0.239 

0.235 


In order to examine the nature of various physical quantities throughout the distribution, we have 
considered a particular star 4U 1820-30 for which mass M = 1.58M 0 , radius R = 9.1 km, the 
physical parameter po — 1.08 and the geometric parameter k = 3.1. We have shown the variation 
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of density and pressures in Figure 2 and Figure 3, respectively. It is observed that the transverse 
pressure p± is less than the radial pressure for r in the range 0 < r < 2.78903. Subsequently p ± 
dominates p r in the region 2.78903 < r < 9.1. The radial pressure p r vanishes at r — 9. 1. In Figure 
4, we have shown the variation of anisotropy throughout the distribution. The variations of sound 
speed in the radial and transverse directions are shown in Figure 5. From Figure 6 , it is evident 
that the strong energy condition, p — p r — 2p ± > 0 is satisfied throughout the distribution. Though 
we have not assumed any explicit expression for the EOS in our model, we have shown the nature 
of variation of pressures p r and p± against density in Figure 7. For a relativistic model to be 
stable in its region of validity, we must have the adiabatic index r>|. The variation F against 
radius is shown in Figure 8 . It is clear from Figure 8 that r > 3 throughout the star. The variation 
of gravitational red shift, z(r) — \Je ~ v M — 1 in the radial direction is shown in Figure 9. It is 
easy to note that the red shift is monotonically decreasing function from the centre to boundary. 
Further, the red shift at the centre z c and on the boundary zr are both positive and finite. 



r(Km) 


Fig. 2.— Density Ys Radius 
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Fig. 3.— Radial and Transverse Pressures Vs Radius 

5. Discussion 

Spherical distribution of matter on pseudo-spheroidal spacetimes have been studied by 
a number of researchers in the recent past Tikekar and Thomas (1998); Tikekar (1999); 
Tikekar and Thomas (2005); Thomas et al. (2005); Thomas and Ratanpal (2007); Paul el al. 

(201 1); Chattopadhyay and Paul (2010); Chattopadhyay et al. (2012). In this paper, we have 
obtained a new class of solutions to Einstein’s field equations for a spherically symmetric 
anisotropic distribution of matter and have shown that our model can fit to the observational data 
of a number of well studied pulsars ( Gangopadhyay et al. (2013)). On assuming a particular form 
of radial pressure and on the basis of elementary criteria for physical acceptability of a compact 
spherically symmetric distribution of matter, we have obtained the bounds for the physical as well 
as geometric parameters of the model. It is found that our model can accommodate a number of 
pulsars like 4U 1820-30, PSR J1903+327, 4U 1608-52, Vela X-l, PSR J1614-2230, SMC X-4 
and Cen X-3. We also have studied, in detail, a particular pulsar 4U 1820-30, and have shown 
graphically the profile of different physical quantities throughout the distribution. In short, study 
of compact stars on the background of pseudo-spheroidal spacetime is highly interesting in the 
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0 2 4 6 8 


r (Km) 

Fig. 4.— Anisotropy Vs Radius 

sense that it generates models compatible with observational data and at the same time having a 
definite 3-space geometry, namely, pseudo-spheroidal geometry which many other spacetimes 
may not possess. 
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Fig. 5.— (Sound Speed) 2 Vs Radius 



Fig. 6.— Strong Energy Condition Vs Radius 
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Fig. 7.— Equation of State 



r (Km) 


Fig. 8.— Adiabatic Index Vs Radius 
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r (Km) 


Fig. 9.— Gravitational Red Shift 
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